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The purpose of this monograph is to organize and present, for effective use in spacecraft develop- 
ment, the significant experience and knowledge accumulated in development and operational 
programs to date. It reviews and assesses current design practices, and from them establishes firm 
guidance for achieving greater consistency in design, increased reliability in the end product, and 
greater efficiency in the design effort. The monograph is organized into three major sections that 
are preceded by a brief introduction and complemented by a set of references* 

The State of the Art, section 2, reviews and discusses the total design problem, and identifies 
which design elements are involved in successful design. It describes succinctly the current tech- 
nology pertaining to these elements. When detailed information is required, the best available 
references are cited. This section serves as a survey of the subject that provides background 
material and prepares a proper technological base for the Design Criteria and Recommended 
Practices. 

The Design Criteria , shown in section 3, state clearly and briefly what rule, guide, limitation, or 
standard must be considered for each essential design element to insure successful design. The 
Design Criteria can serve effectively as a checklist for the project manager to use in guiding a 
design or in assessing its adequacy. 

The Recommended Practices , shown in section 4, state how to satisfy each of the criteria. When- 
ever possible, the best procedure is described; when this cannot be done concisely, appropriate 
references are provided. The Recommended Practices , in conjunction with the Design Criteria , 
indicate how successful design may be achieved. 

The design criteria monograph is not intended to be a design handbook, a set of specifications, or 
a design manual. It is a summary and a systematic ordering of the large and loosely organized 
body of existing successful design techniques and practices. Its value and its merit should be 
judged on how effectively it makes that material available and useful to the user. 
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FOREWORD 


NASA experience has indicated a need for uniform design criteria for space vehicles. Accord- 
ingly, criteria are being developed in the following areas of technology: 

Environment 

Structures 

Guidance and Control 
Chemical Propulsion 

Individual components of this work will be issued as separate monographs as soon as they are 
completed. This document, “Spacecraft Aerodynamic Torques,” is one such monograph. A list of 
all monographs in this series issued prior to this one can be found on the last page of this document. 

These monographs are to be regarded as guides to design and not as NASA requirements, except 
as may be specified in formal project specifications. It is expected, however, that the criteria 
sections of these documents, revised as experience may indicate to be desirable, eventually will 
be uniformly applied to the design of NASA space vehicles. 

This monograph was ’prepared under the cognizance of the NASA Electronics Research Center 
and published by the Jet Propulsion Laboratory. Robert Lyle and Pericles Stabekis of Exotech, 
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University of Illinois, and R. Passamaneck of the Jet Propulsion Laboratory. 
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Contributions in the form of design and development practices were provided by many other j 

engineers of NASA and the aerospace industry. I 

Comments concerning the technical content of these monographs will be welcomed by the Na- 
tional Aeronautics and Space Administration, Office of Advanced Research and Technology (Code 
RE), Washington, D.C. 20546. ; 

January 1971 J 
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SPACECRAFT 

AERODYNAMIC TORQUES 


1 . INTRODUCTION 

In the design of spacecraft attitude control systems all torques which tend to disturb the attitude 
of a spacecraft must be considered. One of these is the aerodynamic torque resulting from the 
impingement of atmospheric gas molecules on spacecraft surfaces. 

Determination of the atmospheric torque requires knowledge regarding 

(1) The atmosphere and its motion resulting from the Earth’s rotation 

(2) The spacecraft aerodynamic and mass characteristics 

(3) The interaction of the spacecraft with the atmosphere and the relative velocity of the 
spacecraft with the atmosphere 

Aerodynamic torques diminish rapidly as the orbital altitude is increased. In general, for space- 
craft in Earth orbit, the radiation force on a given surface becomes more significant than the 
aerodynamic force at orbital altitudes above 1000 km. Conversely, below 600 km aerodynamic 
forces usually predominate. Between 600 and 1000 km, torques resulting from solar radiation 
and aerodynamic forces are likely to be of the same order of magnitude. However, the aero- 
dynamic forces on the spacecraft, from which the aerodynamic torques originate, are the main 
factors in the determination of spacecraft orbital lifetime. 

The aerodynamic torque may be an important consideration in the determination of spacecraft 
attitude motion, control actuator sizing, and fuel requirements. Control or minimization of atmo- 
spheric drag torques requires that attention be given to the geometry and characteristics of the 
external surface of the spacecraft, particularly the configuration and surface properties of extended 
structures such as large solar arrays, antennas, and booms. 

The scope of this monograph is limited to the assessment and accommodation of disturbance 
torques due to the interaction of spacecraft with the atmosphere in comparitively long duration 
orbits. This typically means flight in the free molecular flow regime which occurs at Earth alti- 
tudes of 120 km or greater. Atmospheric models for Earth, Mars, and Venus recommended for 
use are given in other monographs (refs. 1 to 4). 

Criteria and recommended practices for the assessment and control of magnetic, gravitational, 
radiation, and mass expulsion disturbance torques are given in other monographs (refs. 4, 5, 6, and 
7, respectively). The relative magnitudes of environmental disturbance torques on an Earth 
satellite are briefly presented in appendix A. Dynamic effects resulting from the coupling and 
interaction of disturbance torque effects with the spacecraft control system are treated in 
reference 8. 
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2. STATE OF THE ART 
2.1 General 

The interaction between a body and an atmosphere through which it is moving has been studied 
for centuries and the accumulated knowledge composes the technology of aerodynamics. The 
particular part of this technology that has importance to problems of attitude stabilization of 
spacecraft relates to the disturbance torque which may be caused by the interaction of a vehicle 
with the rarefied atmosphere at orbital altitudes. 

At orbital altitudes, the interaction of the vehicle and the very rarefied atmosphere can best be 
characterized by the free-molecular flow regime of gas dynamics, i.e., when the molecular mean 
free path is much greater than a characteristic spacecraft dimension. For spacecraft of approxi- 
mately 1-m cross-sectional diameter, this condition is satisfied for altitudes above 120 km. 


2.2 Historical Background 

Rarefied gas dynamics has been the subject of numerous investigations since the time of Maxwell. 
The early studies were limited to cases of very low speeds and, generally, to "internal” flow geom- 
etries associated typically with vacuum installations. Since about 1946, the possibility of flight 
at high altitudes and high speeds prompted renewed interest and activity in the field. A number 
of experimental and theoretical results were obtained, with emphasis on the aerodynamic problem 
of high-speed flow of a rarefied gas past a body (ref. 9). 

In rarefied gas flow, the length of the molecular mean free path, A, is comparable to some signifi- 
cant flow field dimension L For the considerations of this monograph, l is a characteristic dimen- 
sion of the spacecraft. The dimensionless ratio A / i is called the Knudsen number and is denoted 
by K. The Knudsen number serves as a criterion for the division into various flow regimes. Thus, for 
K< 10~ 2 the flow is classified as "continuum flow;” for 10 2 < K<10 _1 , "slip flow;” for 10 _1 < K<10, 
"transition flow;” and for K> 10, "free-molecular flow.” In some cases, primarily blunt bodies, free- 
molecular flow may exist down to K = 3. 

Early theories concerning interactions of impinging particles on surfaces have been based either 
on a quantum or classical formulation of the interaction. Representative of the quantum approach 
are the works of Jackson (ref. 10) and Devonshire (ref. 11). These earlier wave mechanics theories, 
however, suffer from the inability to deal successfully with realistic interactions, particularly of 
the heavier gases. This was pointed out by Zwanzig (ref. 12) and Goodman (ref. 13), whose works 
are representative of the classical treatment of the problem. 

Observations of rarefied gas flows relating to problems of molecular interaction at the surface 
were first made near the turn of the century. 
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In early experiments by Millikan (ref. 14) the drag on a rotating cylinder due to the impact of 
molecules of rarefied gas was treated in terms of a parameter representing the exchange of tan- 
gental momentum between the cylinder wall and the incident particles. This parameter was related 
to the Maxwellian concepts of specular and diffuse reflection of light and is identified as the 
coefficient of tangental momentum exchange <r t . 

Another parameter introduced in this period was the accommodation coefficient for thermal energy 
a. This coefficient was introduced by Knudsen in his efforts to interpret experiments on the transfer 
of thermal energy from a fine wire (also discussed in ref. 14). Both Millikan’s and Knudsen’s 
experiments were performed under large Knudsen number conditions (free-molecular flow). The 
use of the two coefficients, a t and a, proved so successful that all free-molecule kinetic theory 
involving surface phenomena remains formulated in these terms. A third coefficient, <r„, the coeffi- 
cient of normal momentum exchange, was first introduced by Schaaf and Chambre (ref. 15) in 
1958. These terms <r f , a„, and a are defined in appendix B. 


2.3 Flight Experience 

The early “paddlewheel” satellites (with solar cell arrays deployed like large windmill airfoils), 
produced some of the first flight experiences with aerodynamic disturbance torques. These satel- 
lites were designed to study the magnetosphere and the solar wind, and it was only coincidental 
that they revealed information about the interaction of the satellite’s surfaces with the atmosphere. 

When the first of these satellites. Explorer 6, was designed, it was expected that the momentum 
of the atmospheric molecules reemitted by its surface would be described by Maxwell’s classical 
model, with approximately 95 percent of the molecules diffusely reflected. A postflight analysis, 
however, revealed that the spin had decayed three times faster than expected, and one order of 
magnitude faster than would have occurred if the reemission had been completely diffused and 
accommodated (ref. 16). The Explorer 6 results were not published and, as a result, designers 
continued to underestimate the aerodynamic torques on paddlewheel satellites. This lack of docu- 
mentation resulted in similar problems for Ariel 2, many of whose experiments became inoperative 
after 3 months because of the rapid decay of the spin, and the consequent attitude drift (ref. 17). 
After the experience with Ariel 2, designers used higher perigee altitudes for paddlewheel satel- 
lites, thus reducing the aerodynamic torques. 

Explorer 17 also experienced a decay of its spin rate attributed to the propeller mode aerodynamic 
effect. Because the despin torque was due entirely to reemission momentum and the orbit decay 
was caused by the combination of incident and reemitted momenta, the data from the combined 
decay effects has been used to study surface interaction effects independent of the actual atmo- 
spheric density (ref. 17). 

The Orbiting Solar Observatory satellite (OSO 1) was subjected to aerodynamic forces which acted 
to roll the spacecraft about the solar vector and to precess the spin axis away from the desired 
pitch position, i.e., normal to the solar vector (ref. 18). This effect was corrected by the attitude 
control system in flight operations. 
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Little is known about the Soviet flight experience; however, in commenting upon the spin decay 
of Sputnik 1, Warwick attributed the reduction of spin velocity to atmospheric friction acting on 
the long antennas of the satellite (ref. 19, p. 155). 


2.4 Evaluation of Surface Interaction 
With the Atmosphere 

The aerodynamic disturbance torque acting on an orbiting spacecraft is the net sum of the indi- 
vidual torques acting about its center of mass. In many cases, the aerodynamic torque is insignifi- 
cant when compared to other disturbance torques. For this reason, simple calculations are often 
made to obtain conservative estimates of the aerodynamic torque levels. Should these simple 
calculations indicate that the net aerodynamic torque is significant when compared to other dis- 
turbance torques, then more complex, but more accurate, methods are employed. 


2.4.1 Aerodynamic Force 

In the rarefied atmosphere at orbital altitudes, the gas molecules that hit the spacecraft are 
reemitted and travel far before colliding with other molecules. In this regime of free molecular 
flow gas dynamics, the effect of the reemitted particles on the incident stream can be neglected, 
at least so far as subsequent effects on the spacecraft are concerned. Therefore, for the purposes 
of this monograph, the incident flow is considered to be undisturbed by the presence of the 
spacecraft. 

This noninteraction of the incident and reemitted particles allow the net aerodynamic torque to 
be calculated by summing up the torque contributions of each of the spacecraft elements. Thus, 
the vehicle may be dissected into simple subshapes to facilitate estimation of the aerodynamic 
forces. Summing these forces gives the result for the entire spacecraft. 

However, before this dissection process is undertaken, an approximation of the total aerodynamic 
force acting on the vehicle is often used to obtain a conservative estimate of the aerodynamic 
torque. The gross value of the aerodynamic force for each applicable vehicle orientation may be 
obtained using the following expression: 

F = Ic„pV 2 A (2-1) 

where 

F = total aerodynamic force 
p = atmospheric density 

A — projected area of spacecraft element normal to the incident flow 
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V = spacecraft velocity 

C D = drag coefBcient, which for the purposes of a conservative estimate of the force is gen- 
erally assumed to have a value of 2.6 

It should be noted that for the calculation of spacecraft lifetime a more accurate drag coefficient 
must be used. For calculations of spacecraft lifetimes, see reference 20. 

A better approximation of the total aerodynamic force can be obtained by using simplified 
particle/surface interaction models, based on free-molecular flow theory. Applying one of the most 
commonly used interaction models (ref. 21), the aerodynamic force on an elemental area is: 

dF - ~ pV 2 [(2 - o- n - a t ) (e„ *e„) 2 e n + <r f (e P *e„) e l? ] dA (2-2) 


where 

o-„ = normal momentum exchange coefficient 
<r f = tangential momentum exchange coefficient 
e v = unit spacecraft velocity vector 
e n = outward unit vector normal to dA 

The linear relationship between dF and the atmospheric density p indicated by the above expres- 
sion, as well as other interaction models, calls for an accurate model of atmospheric density. The 
momentum exchange coefficients are generally considered to be functions of the surface material 
of the spacecraft. However, for this approximation, an empirical value of 0.8 has been used for 
o t and cr n . 

When the use of such approximations indicates that the aerodynamic force is significant, as 
compared to other environmental forces, a more detailed calculation may be necessary. Such a 
calculation goes beyond empirical estimates of the surface interaction and requires an under- 
standing of the basic physics of particle interactions that depend on surface roughness, tempera- 
ture, and composition of both the gas and the surface. The differential force dF for an arbitrarily 
oriented surface element dA, in terms of direction cosines (ref. 22), is: 

dF = j-^pV 2 ^«r,(efc + v t) + (2 - <r„) £p][f(l + erf fS) + -^=- exp (-£ 2 S 2 )] 

+ p(l + erf £S) + -^[nr (1 + erf(S) + ■F ex P(-^ S2 )])} dA 

(2-3) 

where 


dF ~ component of the force in the direction specified by the direction cosines k , p, t 

c, f , r} = direction cosines specifying the orientation of the incident velocity vector V 
(see fig. 1) 
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k } p,t = direction cosines specifying the direction of dF (see fig. 2) 


S = 


molecular speed ratio = 



R = gas constant 
m = molecular weight 
Ti = temperature of incident molecules 
T w — surface temperature 


2 fz s 

erf (fS) = — — ( e“<*> a dx 

J o 


exp(-£ 2 S 2 ) = 


The direction cosines fc, p, t for the aerodynamic force are the cosines between the force vector 
and the x, y, z coordinates, respectively. Similarly, the direction cosines e, £, rj for the velocity 
vector are the cosines between the velocity vector and the x, y, z coordinates, respectively. The 
reference frame, the x, y, z coordinate system, is arbitrary and need not be orthogonal. The 
convenience of using an orthogonal coordinate system is that the three components of velocity or 
force (if desired) are mutually perpendicular and independent of each other. The orientation 
of the coordinate system on the spacecraft is also arbitrary; however, if the system is oriented so 
that the f/-axis is normal to the body surface at a convenient point, then the angle of attack, 0, is 
the angle whose cosine is equal to f . 

The total force acting on the spacecraft can be obtained by integrating dF over the entire space- 
craft surface. It should be noted that dF is the net force acting on dA , combining the force com- 
ponents due to both incident and reemitted molecules. Equation (2-3) shows that the net force 
generally depends on the normal momentum exchange coefficient <r n as well as on the ratio of 
surface temperature T w to the incident temperature Ti and hence on the thermal accommoda- 
tion coefficient a. In practice, it is frequently assumed that the reflection is completely diffuse 
(<t„ — 1) and that the surface temperature T w is uniform. 

In a more realistic calculation, the surface temperature distribution due to effects of accommoda- 
tion and reemission should be determined taking into account the internal conduction and external 
radiation characteristics of the body. This is done by establishing a heat balance for the spacecraft 
which will give the net heat flow rate. For a complete discussion of the heat flow due to accom- 
modation and reemission as well as how to apply it, see ref. 9. The force should then be inte- 
grated over the surface using appropriate local values of the momentum exchange coefficients and 
surface temperature. 
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2.4.2 Aerodynamic Torque 

When the expression given by equation (2-1) is used to evaluate the aerodynamic force, the aero- 
dynamic torque is, in turn, estimated by the following expression: 


L a = IF 


(2-4) 


where 


L a — total aerodynamic torque 
i = moment arm 


and 


F is as given by equation ( 2-1 ) 

For the purposes of conservatively estimating the torque from this model, the moment arm i is 
assumed to be at least one-third of the spacecraft's maximum dimension, including all appendages, 
even if the spacecraft is symmetrical, and greater if the projected area centroid and center of mass 
are farther apart. The center of pressure is a useful concept that may be applied to estimate 
aerodynamic torques at the outset of a problem. The center of projected area may be used as the 
line of action of the aerodynamic force as a first approximation, recognizing the variable way in 
which the flow interacts with the surface at different locations. A conservative estimate for the 
moment arm is used (i.e., no less than one-third the vehicle's maximum dimension). 

Evaluation of the aerodynamic torque through the application of the simplified particle-surface 
interaction model defined by equation (2-2) or the more exact expression for the aerodynamic 
force given by equation (2-3) involves the calculation of the torque contribution from the various 
surfaces of the spacecraft. The usual procedure is to approximate these surfaces by means of 
simple geometric shapes (planes, cylinders, cones, spheres, etc.) so that the resulting surface inte- 
grals may be readily evaluated. If Jb is the vector from the spacecraft center of mass to an 
infinitesimal area dA { on the surface acted upon by a differential force then the torque con- 
tribution from this surface is 


L = 1* X dFi (2-5) 

J surface 

The total aerodynamic torque on the spacecraft is then obtained from the vector sum of the 
torques on the elementary shapes which approximate the spacecraft, i.e., 

L„ - £ L,,; (2-6) 

i -1 

where n is the number of surfaces chosen to model the spacecraft. 

The torque calculation can be complicated if successive interaction of the flow from one surface 
onto another is considered. This may be an important consideration for low Knudsen numbers and 
concave surfaces. 
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Torques on paddlewheel satellites have shown the importance of considering more accurate 
models of the aerodynamic interaction than the noninteracting model for incident and reemitted 
particles. For crude estimates, a process similar to the modelings used with equations (2-1) and 
(2-4) may be used. The force term can be taken as (1 4- r) times the force of the incident mo- 
mentum change, where r is the ratio of the speed of reemitted molecules to the speed of incident 
molecules. The torque is then evaluated assuming an appropriate moment arm. This torque is 
contained as a component of L a in equation (2-6); therefore, equation (2-6) is valid when the more 
exact models defined by equations (2-2) and (2-3) are employed. 


2.4.3 Aerodynamic Coefficients in Free Molecular Flow 

It is generally convenient to express the results of aerodynamic force and torque analyses in terms 
of a drag coefficient, C D , 1 and a moment coefficient, C M . These coefficients are given by 

= ^ (2-7) 

T pV->A 

and 


c u = — r — — (2-8) 

JPV'-AI 

Since a body in free-molecular flow does not disturb the flow, a complicated shape can be 
resolved into simple parts, and the contributions of each of these parts can be added together to 
obtain the coefficients for the entire spacecraft. Thus, it is generally sufficient to determine the 
aerodynamic coefficients for simple shapes, i.e., flat plates, cylinders, cones, etc. 

When the speed of a spacecraft is large, relative to the random thermal motion of the molecules, 
the flow may be described as hyperthermal. Care must be exercised when hyperthermal flow is 
assumed and the random thermal motion of the molecules neglected for calculating aerodynamic 
coefficients for approximate models. Unacceptable errors in drag and moment coefficients will 
result when the spacecraft has large surfaces that are almost parallel to the flow. These errors arise 
from the tangential momentum exchange which is due to the random thermal motion of the 
molecules. The error introduced by this assumption is discussed further in references 23 and 24. 
Table 1 (ref. 23) lists the drag coefficients of spacecraft in hyperthermal free-molecular flow for 
certain simple shapes. The drag coefficients are listed in terms of the ratio r of the speed of 
reemitted molecules v r to the speed of incident molecules v*. The dependence of the drag coeffi- 
cient on the thermal accommodation coefficient a is illustrated in figure 3 for a flat plate in 
hyperthermal free-molecular flow with the angle of attack 0 as the parameter. 


1 The drag coefficient is often expressed in terms of an axial force coefficient C 4 and a normal force coefficient C y : 
C d~C a cos 0 + C v sin 6, where 6 is the spacecraft angle of attack with respect to V. 
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TABLE 1. Drag Coefficient in Hyperthermal Free-Molecule Flow (ref. 23) 


Shape 

Drag coefficient, C D (based on projected area 
perpendicular to direction of motion ) 


Diffuse reemission 

Specular reflection 

Flat plate ( normal to flow ) 

2 ( i+ f') 

2(1 +r) 

Flat plate ( at angle of attack 0 ) 

2 ^1 + r sin 0^ 

2 (1 — r cos 20) 

Sphere 

2 ( 1 + 1 ') 

2 

Cylinder ( perpendicular to flow ) 

2 ( 1 + i') 

*( 1 + T ') 

Cone of semi vert ex angle $ with vertex 
forwards and axis parallel to flow 

2 (l + -|- rsin 

2(1 — r cos 2^ ) 


When a more exact calculation of the aerodynamic coefficients is needed, the random thermal 
motion of the atmospheric molecules should be included. Reference 22 gives a detailed analysis 
of such treatment and presents aerodynamic coefficients obtained for certain simple shapes. These 
results also include the effects of self -shading, such as the shading of the back side of a cylinder 
by its front side. Self-shading effects are also included in equation (2-3), but not the shading of 
one subshape by another. 

The shading mentioned above is the part of the spacecraft that has been blocked from the incident 
flow due to another part of the spacecraft. There are two commonly used methods for estimating 
shading. One method is to use “light-ray” shading. In this method, the “light” is assumed to be a 
beam parallel to the flow velocity of the incident molecules. Any subshape or part thereof 
that falls in the shadow cast by another subshape is assumed completely shadowed from the 
incident flow and, hence, yields no contribution to the drag or moment coefficients. The second 
method considers the random motion of the molecules superimposed upon the flow velocity. Con- 
sequently, the molecules penetrate into the “light-ray” shadow as described above. The amount 
of this penetration is approximated by the angle tan* 1 1/S, as shown in figure 4. As S— » oo, the 
penetration is nonexistent and ‘light-ray” shading results, thus light-ray shading is valid for a very 
large speed ratio. As S-» 1, the penetration angle approaches 45° and should be considered as 
the lower limit of validity since filling of the shaded region can occur through the random motion 
of the molecules. 
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2.5 The Atmosphere at Orbital Altitudes 

Atmospheric conditions encountered by a spacecraft in an Earth orbit are of great importance 
to the study of aerodynamic forces and torques. The torques resulting from the interaction be- 
tween the spacecraft and the atmosphere are determined to a large extent by atmospheric density. 
It is essential that accurate atmospheric models be used which provide atmospheric density as well 
as chemical composition and temperature. Two such models for the Earths atmosphere are pre- 
sented in reference 1 — the Jacchia’s modified model, developed in 1967, and the Quick-Look 
Density model The former is a detailed computerized model which can describe the atmosphere 
from altitudes of 120 to 1000 km for any time and spatial location. Figure 5 shows four atmospheric 
density profiles provided by this model and also shows the importance of the sunspot cycle 
on density. The Quick-Look Density model was developed mainly for preliminary mission planning 
and provides an estimate of the expected atmospheric density at the planned orbital altitude 
without the use of a computer. 


2.6 Vehicle Characteristics 


The interaction of a spacecraft with the atmosphere is largely affected by several of the space- 
craft’s characteristics. Thus, there is a functional dependence of the aerodynamic torque on the 
vehicle’s configuration, the attitude angles, the location of the center of mass relative to the center 
of pressure, and the surface properties. 


The shape of the spacecraft and the attitude angles, which relate the orientation of the geometric 
axes to the flow vector, are necessary for torque calculations since the angle of attack required in 
the expressions given for the aerodynamic force and torque is directly related to them. Typically, 
the nominal attitude angles (assuming perfect control) are selected. 


The location of the spacecraft center of mass relative to the center of pressure may be useful when 
considering aerodynamic torques. The center of pressure can be calculated from the equation 


X cp 


X 


cm 



(2-9) 


where 


X cp — the distance of the center of pressure from a designated reference point 

X cm = the distance of the center of mass from the same reference point 

La x = the aerodynamic torque about the x-axis 

F„ x = the aerodynamic force component normal to the flow; i.e., 

[(XF y Y + (2F*) 2 p 
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C M = the moment coefficient 

C N — the normal force coefficient 

JL ref = the reference length (moment arm) 

In equation (2-9) it should be noted that X cp , X cm , and the reference point must be colinear. To 
avoid confusion, the reference point may be taken as X cm . Equation (2-9) shows that, to calculate 
X cp , it is first necessary to calculate L ax and Fn x ■ Thus, the center of pressure is a derived result 
and, generally, is not calculated in practice. It is evident from this equation that it cannot repre- 
sent torque components parallel to the incident force. Therefore if needed, Y c „ — Y cm — L aj /F„ v 
and Tjcp 2 cm = L a JF nz can be calculated, where L ay , L Qz , F n , y , and F Uz have analogous defini- 
tions to L ax and F Hx as given above. 

Other spacecraft surface characteristics of importance include the atomic composition of the 
outermost layer as well as its temperature distribution. The thermal accommodation coefficient a 
is dependent on the ratio of the mass of the incident gas atoms to the mass of the surface atoms. 
Thus, a rigorous analysis of the interaction is dependent on knowledge of the surface layer of 
atoms. For a critical survey of thermal accommodation coefficients, see reference 25. 


2.7 Testing 

There are no presently available testing techniques for determining, prior to flight, the effects 
of spacecraft interaction with the environment in which it will operate. Particular areas which 
have hindered the development of adequate testing techniques include simulation of the environ- 
ment at orbital altitude and measurement of the distribution of reflected molecule flux and energy. 

Experimental research of molecule/surface interactions has been done with nozzle-generated 
molecular beams to produce a high-speed ratio collisionless beam (ref. 26). Efforts to increase the 
stagnation temperature and lower the mean molecular weight of the source gas mixture, so that 
Earth orbital velocity beams can be produced, are continuing. Molecular beam experimentation, 
performed in the late 1950 s and early 1960 > s with apparatus not quite capable of simulating the 
orbital environment, yielded values for the accommodation coefficients lower than those indicated 
by past flight experience (ref. 27). 

Recent efforts for better Earth orbital simulation include the development of a partially dissociated 
atomic oxygen beam (ref. 26). The development of this beam represents a major step toward 
complete compositional simulation since atomic oxygen is the predominant constituent of the 
atmosphere for the orbital altitudes where aerodynamic effects are strong. 
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2.8 Summary 


Flight experience has shown that aerodynamic torques can significantly disturb the attitude of a 
spacecraft, particularly by affecting the spin rate, and therefore must be considered in the design 
of the spacecrafts attitude control system. Aerodynamic torque is an important consideration in 
the design of large spacecraft and gravity-stabilized satellites with a perigee altitude ranging 
between 120 to 1000 km. Aerodynamic forces in this altitude range are calculated using free- 
molecular flow theory. At altitudes above 1000 km the radiative force dominates compared with 
aerodynamic torques. Down to altitudes of 600 km, the radiative force is of the same order of 
magnitude as the aerodynamic force and, therefore, must be included in the calculations. Below 
an altitude of 120 km, the aerodynamic forces become large and, hence, the lifetime of the space- 
craft would be very short. Also, below 120 km, the assumption of free-molecular flow theory may 
not be valid. Analytical techniques for determining the aerodynamic torque require several simpli- 
fying assumptions and approximations to make the problem solvable by numerical techniques. The 
results, however, are sufficiently accurate to enable prediction of a spacecraft’s performance and 
lifetime. 


3. DESIGN CRITERIA 

Disturbance torques arising from aerodynamic forces acting on spacecraft surfaces must be 
accounted for in the design of attitude control systems. It should be demonstrated that the aero- 
dynamic torques acting in combination with all other disturbance torques do not degrade the 
performance of the attitude control system. When it is determined that aerodynamic disturbance 
torques may be an important factor in the attitude control system performance, procedures for 
minimizing or accommodating associated torques should be initiated and followed during the 
spacecraft design, development, and fabrication. 


3.1 Aerodynamic Torque Analysis 

The aerodynamic torque should be determined at sufficient points in the orbit and times of the 
year to generate a profile of the torque under all operating conditions. Care should be taken not 
to eliminate calculation of maximum and minimum torque values that would occur at apogee 
and perigee or at given times of the year. 
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3.1.1 Environment 


The atmosphere chosen for making calculations of aerodynamic torques on orbiting spacecraft 
should be sufficiently accurate so as not to cause errors greater than the maximum error deter- 
mined by attitude control system performance. A secondary criterion that should be considered 
is the ease with which the atmospheric model may be applied. 


3.1.2 Surface Orientation 

Each part of the spacecraft should be represented by the simplest possible subshape for aero- 
dynamic torque calculations. Only surfaces exposed to the flow (e„ • e n > 0) should be considered 
at any particular time. The flow vector must always be kept in proper orientation to the subshape 
of interest. 


3.1.3 Center of Pressure 

The distance between the spacecraft’s center of mass and center of pressure should be estimated 
for various spacecraft orientations relative to the flow vector, since the center of pressure will 
move depending on the orientation. The estimated center of pressure should always be made 
conservative in the light of aerodynamic torque; i.e., the distance between the center of mass 
and center of pressure should be conservatively large. 


3.1.4 Aerodynamic Force 

The net aerodynamic force acting at the spacecraft center of pressure should be determined for 
all expected flight configurations, orientations, and environmental conditions using methods 
described in section 2.4.1. 


3.1.5 Torque Variations 

The variations in aerodynamic torque resulting from shifts in the spacecraft center of mass (caused 
by mass expulsion, deployment of appendages, equipment jettison, etc.), and shifts in the center 
of pressure (caused by relative motions of appendages, changes in surface characteristics, shading, 
etc.) should be evaluated and accounted for to insure proper spacecraft performance. 
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3.2 Evaluation of Aerodynamic Torque Effects 

The evaluation of the effects of aerodynamic torques on spacecraft should include, but not be 
limited to, the following considerations: 

(1) Attitude control system actuator requirements; viz., peak torque, momentum storage, mo- 
mentum transfer, and control system propellant requirements 

(2) Attitude control accuracy 

(3) Structural deflections causing shifts in the center of pressure relative to the center of mass. 

3.3 Control of Aerodynamic Torque 

Procedures should be instituted for the determination and control of surface properties, center of 
mass, and center of pressure whenever aerodynamic torques are dominant or contribute signifi- 
cantly along with other torques to spacecraft attitude disturbances which are large with respect 
to the control system capacity. These procedures should be initiated in the early design phase 
and maintained throughout the development program. 


4. RECOMMENDED PRACTICES 

It is recommended that aerodynamic torques be assessed in the early design phase of spacecraft 
development. Experimental or test techniques to ascertain the magnitude of aerodynamic torques 
on a spacecraft prior to flight do not exist. Therefore, current practices for estimating these 
torques should be based on knowledge accumulated from previous flight experience and on 
approximate analytical techniques. 


4.1 Aerodynamic Torque Analysis 

4.1.1 General Procedure 

The objectives of a torque analysis in the preliminary phase of spacecraft design are (1) to achieve 
a reasonable approximation of the magnitude of the torque, (2) to identify the spacecraft con- 
figurations having the greatest potential for causing attitude stabilization problems, and (3) to 
determine the constraints imposed on the attitude control system for each proposed spacecraft 
configuration. 
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The preliminary analysis should establish guidelines for tradeoffs associated with various design 
configurations to indicate whether more precise analysis is needed. Approximate location of the 
center of mass at this preliminary stage can only be based on gross mass estimates for the main 
spacecraft elements and all major appendages. Similarly, the approximate location of the center 
of pressure can only be based on the probable spacecraft geometry which can be approximated 
by combinations of tetrahedrons, cylinders, spheres, cones, and plates along with estimates of 
their surface properties. 


When aerodynamic torques are of the same order of magnitude as other significant torques which 
cause spacecraft attitude disturbance, or when aerodynamic torques are the singular dominant 
torques, a detailed analysis is necessary. This analysis requires close coordination between the 
control systems, aerodynamic, structural design, thermal design, and material specification groups. 
As the design progresses, the participating groups should be aware of all changes in configuration 
and materials so that the effects of these changes may be evaluated. 


Methods of analysis which treat the aerodynamic torque problem with a level of detail that is 
commensurate with the state of the art are discussed in reference 27. These techniques should 
be pursued if the preliminary worst-case analysis shows that aerodynamic torque may cause 
significant disturbance. 


4.1.2 Characterization of Atmospheric Environment 

Two models are recommended for the Earth’s atmosphere: the modified Jacchia model (1967) 
and the Quick-Look Density model, both presented in detail in reference 1. For the atmospheres 
of Mars and Venus, the models described in references 2 and 3, respectively, are recommended. 


4.1.3 Aerodynamic Force 

The most commonly used techniques for determining the aerodynamic force on a spacecraft 
involve the approximation of the surface configuration by using a number of simple geometrical 
shapes. Self -shading, as well as shading of one subshape by another, should be evaluated through 
adequate shading models, such as the ‘light-ray” and 1/S models presented in section 2.4.3. 

The approximate models defined by equations (2-1) and (2-2) are recommended for calculating 
the aerodynamic force in a worst-case analysis. If the results of this analysis indicate that aero- 
dynamic effects arc significant, it is recommended that a detailed calculation be undertaken using 
the more exact model defined by equation (2-3). 

Calculations using the exact model are complex and the assistance of specialists in aerodynamics 
and heat transfer is recommended. When employing this model, surface temperature distribution 
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should be evaluated with the methods described in reference 10, taking into account the internal 
conduction and external radiation characteristics of the body. In this analysis, the recommended 
practice is to assume that spacecraft surfaces expose mainly oxygen atoms to the atmosphere, 
unless these surfaces are painted or gold-plated. Experimental results give evidence that, for most 
materials, the values of the normal momentum exchange coefficient and the thermal accommoda- 
tion coefficient is between 0.8 and 1.0 (ref. 28). These values, however, may be modified in view 
of the on-going efforts for better testing techniques. 


4.1.4 Aerodynamic Torque Acting on the Spacecraft 

The computation of the torque contribution for each surface of the spacecraft, approximated by 
simple geometric shapes, allows calculation of the total torque acting about the spacecraft center 
of mass from the equation 


L = ix ; = £[!,, x FiJ = if XdF.j 

where 

L a ~ aerodynamic torque about the spacecraft center of mass 

L ai = aerodynamic torque about the spacecraft center of mass due to the i th subsurface 

J P i — vector from the spacecraft center of mass to the center of pressure of the i th 
subsurface 

Fi = aerodynamic force acting at the center of pressure on the i th surface 
(expressions for Fi are presented in section 2.4.1) 

When the approximate model defined by equation (2-1) is used for torque calculations, the 
moment arm is taken to be at least one-third of the spacecraft's maximum dimension. This is a 
conservative approximation for the distance between the vehicle s center of mass and the center 
of pressure. The latter becomes a derived result when the more exact model is employed and is 
not a requisite input to the analysis. 


4.1.5 Variation of Aerodynamic Torque 

The dependence of the total aerodynamic torque on various parameters requires that a sensitivity 
analysis be performed to determine the magnitude of the change in torque as a function of the 
departure from the nominal values of these parameters. The torque calculations can be repeated 
assuming worst-case departures from nominal values for such parameters as atmospheric density, 
center of mass, surface characteristics, etc. 
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4.2 Minimization of the Effects of Aerodynamic Torques 


Proper design of the spacecraft can reduce aerodynamic torques to manageable levels. An appli- 
cable technique is to balance the aerodynamic forces about the spacecraft center of mass by 
designing the exposed surfaces so that the proper area distribution, relative orientation of each 
surface, and surface characteristics are achieved. 

For spacecraft employing extensible structures, the orientation that minimizes the aerodynamic 
torque should be established. In general, in an optimum configuration the extensible structures, 
such as solar panels, are coplanar rather than staggered to prevent the possibility of unbal- 
anced torques in an orientation other than the design one. The distinct advantage of choosing 
coplanar panels for attitude control considerations was recently confirmed during the design of 
the Small Astronomy Satellite at the Applied Physics Laboratory, Johns Hopkins University. 

The stabilizing torques for passively stabilized spacecraft are typically smaller than those encoun- 
tered on actively stabilized spacecraft and, therefore, passively stabilized spacecraft are more 
susceptible to disturbances caused by aerodynamic torques. 

A recommended practice for reducing the aerodynamic torque for spacecraft is to reduce the 
eccentricity and increase the perigee altitude, provided that such increase does not affect the 
mission requirements. Orientation of appendages which would allow their torque contributions to 
be additive, such as solar panels arranged in a paddlewheel orientation, should be avoided. 


4.2.1 Spin-Stabilized Spacecraft 

Many spin-stabilized satellites do not have the capability of applying control torques about their 
spin axis, but use inertia to maintain their spin rate. These satellites are quite sensitive to altera- 
tions of their spin rates by aerodynamic torques. If there are surfaces which face the flow over 
part of the orbit, a reduction in spin rate can occur. 


4.2.2 Gravitationally Stabilized Spacecraft 

Gravitationally stabilized satellites orbiting at low altitudes (viz., below 1000 km) are sensitive to 
aerodynamic forces because of the long boom configuration needed to develop the comparatively 
small magnitude gravitational restoring torque. A torque must be added by the control system in 
order to maintain attitude. If not, the vehicle will rotate toward the position in which the gravity 
gradient torque balances the aerodynamic torque. The attitude error that can be tolerated deter- 


20 



mines the minimum altitude of a circular orbit (ref. 29). Eccentric orbits and circular orbits, which 
pass through the diurnal solar bulge in the atmosphere, may have even larger forced response at 
orbital frequency than the steady offset of a constant aerodynamic torque. 


4.2.3 Aerodynamically Stabilized Spacecraft 

When a spacecraft is stabilized by aerodynamic forces (angle of attack 0 = 0) at a given altitude, 
the restraining force is a function of the deviation from 0 = 0 (ref. 29). At altitudes usually con- 
sidered suitable for aerodynamic stabilization, between 120 and 600 km, the aerodynamic forces 
provide negligible damping, so that some means of damping must be provided, such as eddy 
current, viscous fluid damping, or gyros (ref. 30). Whenever feasible, the spacecrafts moments 
of inertia should be designed to be nearly constant about all axes to enhance aerodynamic stabili- 
zation (ref. 29). 


4.2.4 Actively Stabilized Spacecraft 

Aerodynamic torques typically present no difficulty for actively stabilized spacecraft that have 
fuel reserve commensurate with the duration of the mission. To insure that the spacecraft is 
capable of controlling aerodynamic torque disturbances, a calculation of the largest aerodynamic 
torque on the spacecraft should be made. This should be done for the combination of spacecraft 
attitude and orbital altitude that create the greatest disturbance torques. This will insure that the 
actuators have enough authority to overcome the disturbance. A separate calculation should be 
made for the net impulse caused by the disturbance. 

Features can be incorporated into the design of the spacecraft that will operate in low orbits so 
that aerodynamic torques will stabilize rather than disturb the attitude. This can be accomplished 
by adding surfaces to the spacecraft structure so that the optimum distribution of aerodynamic 
forces about the center of mass can be obtained. 
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APPENDIX A 


RELATIVE MAGNITUDES OF SPACECRAFT 
ENVIRONMENTAL TORQUES 


The relative magnitudes of the spacecraft environmental torques on a typical satellite are shown 
in figure A-l (ref. 31). The spacecraft configuration selected for the comparison of these torques 
is a 1.52-m diameter circular cylinder, 9.14 m in length, having the inertia characteristics shown 
in the figure. The spacecraft is assumed to be in a circular orbit. 

The torque due to cosmic dust was evaluated for the cylinder normal to the stream and a 
separation distance of 0.36 m between the center of mass and center of pressure. Evaluation of the 
gravity gradient torque was carried out for an angle of 1° between the cylinder principal axis 
of inertia (the symmetry axis) and the Earth radius vector. The magnetic torque was computed 
assuming a 1-A current in a single loop of wire around the length of the cylinder oriented to 
give the maximum torque. The aerodynamic torque was computed using free-molecular flow 
theory, and the radiation torque was evaluated using a reflectivity of one. 

It should be noted that all of the torques are linear in the parameter shown in the figure, includ- 
ing gravity gradient up to about 15°. Thus, the gravity gradient torque for $ = 10° can be 
obtained by multiplying the values on the appropriate curve of figure A-l by 10; the aerodynamic 
torque for a separation distance c.g. — c.p. = 1.1 m by multiplying the corresponding curve by 3, 
etc. (ref. 31). 
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Maximum torque, Nm (ft— lb) 



Altitude, km (nautical miles) 

Figure A- 1. -Relative magnitudes of the environmental torques on an Earth satellite. 


24 




APPENDIX B 


DEFINITIONS OF THE ACCOMMODATION AND 
MOMENTUM EXCHANGE COEFFICIENTS 


The thermal accommodation coefficient a is a measure of the degree of equilibrium attained 
between the molecule and the surface before the molecule is reemitted. It is defined as (ref. 22) 


a = 


Ej-Er 

Ei - E w 


where 

Ei — the energy carried to unit area of the surface by the incident molecules 

E r — the energy carried away from the unit area by the reflected molecules 

E w — the energy the reflected molecules would carry away from the surface if they were 
reemitted at the temperature of the surface T w 

For complete thermal accommodation E r = E w and, therefore, a = 1. 

The momentum transferred to the surface is specified by momentum exchange coefficients, one 
for the tangential momentum o- t and one for the normal momentum <r n . They are defined as 
(ref. 22) 

_ Ti — T r 
_ 

Ti ~ T W 

_ Pi ~ Pr 
° u Pi - P« 

where 

Ti = the tangential momentum carried to unit area of the surface by the incident molecules 

T r = the tangential momentum carried away from the unit area by the specularly reflected 
molecules 

t w — the tangential momentum which would be carried away from the surface by the 
diffusely reflected molecules if they were in thermal equilibrium with the surface, 
for a = 1. 

Pi — the normal momentum carried to unit area of the surface by the incident molecules 
P r — the normal momentum carried away from the unit area by the reflected molecules 


P w “ the normal momentum which would be carried away from the surface by the diffusely 
reflected molecules if they were in thermal equilibrium with the surface, for a = 1. 


For completely diffuse reflection, r r = t w = 0 and <r t = 1, regardless of thermal accommodation. 
If any specular reflection occurs, v t depends on the degree of thermal accommodation through r r . 
For complete specular reflection and no thermal accommodation, r r = t* and a t = 0- For complete 
diffuse reflection and complete thermal accommodation, <r n = 1, while for complete specular 
reflection and no thermal accommodation, <r n = 0. 


Specular reflection occurs when the angle of reflection is equal to the angle of incidence. Diffuse 
reflection of molecules includes all molecules that are reflected in other than a specular direction. 
Since diffuse reflection is strongly dependent on the reflective surface and its structure, it is 
not symmetrical about the normal to the surface but may have a definite preferred direction 
(refs. 32 and 33). 
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APPENDIX C 
NOMENCLATURE 


A 
C A 
Co 
C m 
Cv 
dA 
d F 
dQ 

Cn 

F 

F 

K 

k,p,t 

L 

l 

7 

m 

R 

r 

S 

Ti 

T w 

V 

Vi 

V r 

a 


Surface area 

Axial force coeflScient 

Drag coefiBcient 

Moment coefficient 

Normal force coefiBcient 

Differential area 

Differential force vector 

Differential heat flow 

Outward unit vector normal to dA 

Unit spacecraft velocity vector 

Aerodynamic force (scalar) 

Aerodynamic force vector 

Knudsen number 

Direction cosines for aerodynamic force 

Spacecraft characteristic dimension 

Moment arm 

Moment arm vector 

Molecular mass 

Gas constant 

speed ratio, — 

V r 

Molecular speed ratio 
Temperature of incident molecules 
Surface temperature 
Spacecraft velocity 
Velocity of incident molecules 
Velocity of reflected molecules 
Thermal accommodation coefficient 
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y Isentropic exponent 
t, £, 7 1 Direction cosines for velocity vector 
0 Angle of attack 
A Molecular mean free path 
p Atmospheric density 
(T n Normal momentum exchange coefficient 
a t Tangential momentum exchange coefficient 
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NASA SPACE VEHICLE DESIGN CRITERIA 
MONOGRAPHS ISSUED TO DATE 




SP-8001 (Structures) 
SP-8002 (Structures) 
SP-8003 (Structures) 
SP-8004 (Structures) 
SP-8005 (Environment) 
SP-8006 (Structures) 
SP-8007 (Structures) 
SP-8008 (Structures) 
SP-8009 (Structures) 
SP-8010 (Environment) 
SP-8011 (Environment) 
SP-8012 (Structures) 
SP-8013 (Environment) 

SP-8014 (Structures) 

SP-8015 (Guidance and 
Control) 

SP-8016 (Guidance and 
Control) 

SP-8017 (Environment) 

SP-8018 (Guidance and 
Control) 

SP-8019 (Structures) 

SP-8020 (Environment) 

SP-8021 (Environment) 

SP-8023 (Environment) 

SP-8024 (Guidance and 
Control) 

SP-8025 (Chemical 
Propulsion) 


Buffeting During Atmospheric Ascent, revised November 1970 

Flight-Loads Measurements During Launch and Exit, December 1964 

Flutter, Buzz, and Divergence, July 1964 

Panel Flutter, May 1965 

Solar Electromagnetic Radiation, June 1965 

Local Steady Aerodynamic Loads During Launch and Exit, May 1965 

Buckling of Thin-Walled Circular Cylinders, revised August 1968 

Prelaunch Ground Wind Loads, November 1965 

Propellant Slosh Loads, August 1968 

Models of Mars Atmosphere (1967), May 1968 

Models of Venus Atmosphere (1968), December 1968 

Natural Vibration Modal Analysis, September 1968 

Meteoroid Environment Model-1969 (Near Earth to Lunar Surface), 
March 1969 

Entry Thermal Protection, August 1968 

Guidance and Navigation for Entry Vehicles, November 1968 

Effects of Structural Flexibility on Spacecraft Control Systems, April 
1969 

Magnetic Fields-Earth and Extraterrestrial, March 1969 
Spacecraft Magnetic Torques, March 1969 

Buckling of Thin-Walled Truncated Cones, September 1968 

Mars Surface Models (1968), May 1969 

Models of Earth s Atmosphere (120 to 1000 km), May 1969 

Lunar Surface Models, May 1969 

Spacecraft Gravitational Torques, May 1969 

Solid Rocket Motor Metal Cases, April 1970 
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SP-8026 (Guidance and 
Control) 

SP-8027 (Guidance and 
Control) 

SP-8028 (Guidance and 
Control) 

SP-8029 (Structures) 

SP-8031 (Structures) 

SP-8032 (Structures) 

SP-8033 ( Guidance and 
Control) 

SP-8034 (Guidance and 
Control) 

SP-8035 (Structures) 

SP-8036 (Guidance and 
Control) 

SP-8037 (Environment) 

SP-8040 (Structures) 

SP-8046 (Structures) 

SP-8047 (Guidance and 
Control) 



Aerodynamic and Rocket-Exhaust Heating During Launch and Ascent, 
May 1969 

Slosh Suppression, May 1969 

Buckling of Thin-Walled Doubly Curved Shells, August 1969 
Spacecraft Earth Horizon Sensors, December 1969 


Spacecraft Mass Expulsion Torques, December 1969 


Wind Loads During Ascent, June 1970 

Effects of Structural Flexibility on Launch Vehicle Control Systems, 
February 1970 

Assessment and Control of Spacecraft Magnetic Fields, September 1970 

Fracture Control of Metallic Pressure Vessels, May 1970 

Landing Impact Attenuation for Non-Surface-Planing Landers, 

April 1970 

Spacecraft Sun Sensors, June 1970 
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